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j/3 ' Abstract 

A theory is presented showing that cloaking of objects from antiplane 
elastic waves can be achieved by elastic pre-stress of a neo-Hookean nonlin- 
ear elastic material. This approach would appear to eliminate the require- 
C/2 . ment of metamaterials with inhomogeneous anisotropic shear moduli and 

_ ^ ' density. Waves in the pre-stressed medium are bent around the cloaked 

C/3 , region by inducing inhomogeneous stress fields via pre-stress. The equa- 

^~»' tion governing antiplane waves in the pre-stressed medium is equivalent to 

the antiplane equation in an unstressed medium with inhomogeneous and 
anisotropic shear modulus and isotropic scalar mass density. Note how- 
ever that these properties are induced naturally by the pre-stress. Since 
the magnitude of pre-stress can be altered at will, this enables objects of 
^ ■ varying size and shape to be cloaked by placing them inside the fluid-filled 

^O ' deformed cavity region. 

^' 

^vq ' Keywords: Cloaking, antiplane waves, pre-stress, neo-Hookean 

rn ' 

2 ■ 1 Introduction 

In recent years there has been a great deal of interest in theoretical and prac- 
tical issues associated with the cloaking of objects when subject to a variety of 
incident wave fields. The original theory was developed in the context of elec- 
tromagnetic waves and used the coordinate transformation principle (Greenleaf 
H ■ et al. 2003, Leonhardt 2006, Pendry et al. 2006, Cummer et al. 2006). Attention 

. - -' then shifted onto the possibility of cloaking in acoustics (see e.g. Cummer and 

Schurig 2007, Chen and Chan 2007, Cummer et al. 2008, Torrent and Sanchez- 
Dehesa 2008, Norris 2008, 2009), surface waves in fluids (Farhat et al. 2008) and 
elastodynamics (Milton et al. 2006, Brun et al. 2009, Norris and Shuvalov 2011). 
It was shown by Milton et al. 2006 that the latter is more difhcult to achieve 
than the other physical applications mentioned due to the lack of invariance of 
Navier's equations under general coordinate transformations, unlike the equa- 
tions governing electromagnetic and acoustic waves. This explains the relative 



smaller number of studies in elastodynamic cloaking. A special case is that of 
flexural waves in thin plates (Far hat et al. 2009). Furthermore invariance of the 
governing equations can be achieved for more a specific class of transformations 
if assumptions are relaxed on the minor symmetries of the elastic modulus ten- 
sor as was described for the in-plane problem by Brun et al. (2009). Norris and 
Shuvalov (2011) exploited this by using Cosserat materials. 

It is well-known that cloaking can be achieved for antiplane elastic waves 
from a cylindrical region (using a cylindrical cloak) in two dimensions (Brun et 
al. 2009). This can immediately be recognized by the duality between antiplane 
waves and acoustics in this dimension. As has been indicated above and as 
will be shown in section [5] of the present paper, in order to achieve cloaking by 
the method of coordinate transformations, the cloak must have an anisotropic, 
inhomogeneous shear modulus and scalar, inhomogeneous density. Indeed, in- 
homogeneity and anisotropy are standard requirements of cloaking materials. 
Furthermore the mappings that require this inhomogeneity are singular in the 
sense that some of these moduli are unbounded on the inner surface of the cloak 
(Greenleaf et al. 2003, Pendry et al. 2006, Cummer et al. 2006). In reality this 
condition must be relaxed so that they are merely required to be very large near 
the inner region of the cloak thus yielding non-perfect cloaking. Construction of 
the cloaking metamaterial itself is non-trivial, although progress has been made 
in this area in some special cases (see e.g. Schurig et al. 2006, Farhat et al. 2008, 
Zhang et al. 2011). 

Here we present an alternative scheme to produce cloaking from antiplane 
elastic waves by the use of nonlinear elastic pre-stress. The nonlinear pre-stress 
generates an anisotropic and inhomogeneous (incremental) shear modulus but 
the density remains homogeneous. Furthermore, these properties arise naturally 
as a result of the pre-stress and so we are not required to construct a special 
metamaterial cloak. Additionally the cloaking is active in the sense that the 
size of the cloaked region can be altered by modifying the pressure inside this 
region. 

We begin with a small cylindrical cavity of radius A located at the origin of 
a soft incompressible nonlinear elastic medium whose constitutive behaviour is 
described by a neo-Hookean strain energy function (Ogden 1997, 2007). With 
reference to figure [H a line source is located at a distance Rq from the origin 
and angle Qq subtended from the X axis. The cavity radius A is chosen so 
that scattering from the line source is weak. This can always be achieved by 
ensuring that KA <^ 1 where K is the wavenumber of the soft rubbery medium. 
We then increase the pressure inside the cavity (by increasing internal fluid 
pressure for example) so that the radius of the cavity region increases to a and 
inhomogeneous stress fields are thus generated within the vicinity of the cavity. 
Note that for a neo-Hookean medium, we require only p-m/ fJ- ~ 3.5 to generate 
a radius that is 20 times larger than its original size and for soft materials we 
have /i ^ 10^ Pa. Furthermore we assume that the fluid inside the cavity is 
inviscid with zero shear modulus so that the appropriate boundary condition 
for the incoming antiplane shear wave is traction free. 

We show that any object placed inside the inviscid fluid in the inflated cavity 



is cloaked in the sense that for a neo-Hookean medium, scattering coefficients 
associated with the inflated cavity (with Ka ~ 0(1)) are exactly the same 
as those associated with the undeformed medium (i.e. scattering is weak). In 
this sense the proposed scheme has connections with reduction of scattering 
by coating regions with special materials as has been discussed in Alu and 
Engheta (2005) although as we stress again, no special materials are required 
here. The effect in the pre-stressed configuration is that waves are bent around 
the large cylindrical cavity by the induced anisotropy and inhomogeneity and 
so an object placed inside this region is not seen by an observer located an 
appreciable distance away. Here we also mention the recent work by Amirkhizi 
et al. (2010) who also discuss the ability to guide waves by continuous changes 
in anisotropy induced by a continuous change of microstructural properties. 

We have non-perfect cloaking in the sense that the original cavity is of finite 
size. This is analogous to the fact that in the classical (linear) cloaking coor- 
dinate transformation, perfect cloaking is not possible because of the singular 
mapping. However here we note that we may always choose this cavity size to 
ensure that scattering is as weak as possible. Furthermore and perhaps most 
importantly, the inhomogeneous moduli and cloaking material are generated 
naturally by the use of nonlinear pre-stress. The continuous change in inhomo- 
geneity induced by pre-stress also means that dispersive effects are not present, 
as they would be by the use of inhomogeneous metamaterials for example. This 
property is beneficial for broadband cloaking. 

We note that nonlinear elastic pre-stress has been increasingly considered 
in recent times in order to control elastic waves, notably in order to tune stop 
and pass bands of composite materials (Parnell 2007, Bertoldi and Boyce 2008, 
Bigoni et al. 2008, Gei et al. 2009). 

In what follows, in section[5]we consider the classical linear coordinate trans- 
formation used in order to achieve cloaking in the antiplane elastic wave con- 
text. In section |3] we introduce the theory of nonlinear elasticity required in 
order to solve the pre-stress problem, and thus determine an equation which 
links the internal pressure inside the cavity to its deformed radius. Section [¥] 
then describes the incremental wave equation which results from considering 
linearizations about the pre-stressed state. This resulting incremental equation 
is the equation which has equivalent form to the standard equation in cloaking 
metamaterials. In section[5]we show how an explicit solution to this incremental 
equation may be obtained thus illustrating that cloaking is theoretically possible 
by the use of pre-stress. We conclude in section |6] indicating possibilities for the 
general elastodynamic case. Some details are located in appendices for ease of 
reading. 



2 Classical coordinate transformation for antiplane 
waves 

Take an unbounded homogeneous elastic material with shear modulus /i and 
density p and introduce a Cartesian coordinate system {X,Y,Z). Let us also 
consider a cylindrical polar coordinate system (i?, 0, Z) with some origin O, 
whose planar variables are related to the Cartesian coordinate system in the 
usual form X = RconQ^Y = i?sin0. Suppose that there is a time-harmonic 
line source, polarized in the Z direction and located at {Rq, 6o), with circular 
frequency uj and amplitude C (which is a force per unit length in the Z direc- 
tion). This generates time- harmonic elastic antiplane waves with displacement 
U = 3fi[W^(X, Y)ez exp(— jojt)] where e^ is a unit basis vector in the Z direction. 
The displacement W is governed by 

Vx • (mVxW^) + P^'^W - ^5{R - Ro)6{e - Go) (2.1) 

Kq 

where Vx indicates the gradient operation in the "untransformed" frame. 

The coordinate transformation theory of cloaking concerns itself with finding 
transformations x to new coordinates x — xOQ which preserve the governing 
equation of motion but result in modifications to the material properties fj, and 
p. This then allows waves to be guided around specific regions of space. The 
transformation can be quite general, the importance being that it possesses some 
singularity (Norris 2008). The classical linear mapping for a cloak for antiplane 
waves (cf. acoustics) expressed in plane cylindrical polar coordinates, takes the 
form (Pendry et al. 2006) 

r^Ri+Rf ^'^^^' Y e = Q, z = Z, (2.2) 

for < i? < i?2 and the identity mapping for all R > R2 for some chosen 
i?i, i?2 G R such that i?2 < ^0 i-e. the line source remains outside the cloaking 
region. The cloaking region is thus defined by r S [i?i, i?2]. 

As is evident from ()2.2|) we use upper and lower case variables for the untrans- 
formed and transformed problems respectively. Therefore under this mapping 
the form of the governing equation (12. ip remains unchanged for R — r > R2, 
i.e. we have 

Vx • (/iVx^w) + pi^^w = ^Sir - Ro)S{e - Go) (2.3) 

whereas for < R < R2, corresponding to the transformed domain Ri < r < 
i?2, the transformed equation takes the form (in transformed cylindrical polar 
coordinates r,6 — Q) 

\d ( , .dw\ ng(r) d'^w ,, , 9 „ ,„ ., 



where 

(2.5) 

Equation (|2.4p is homogeneous since the Une source is located in the outer do- 
main. Note rather importantly that two effects are thus necessary for cloaking 
of this form. Firstly both the shear modulus and density must be inhomoge- 
neous. Secondly the shear modulus must be anisotropic. Finally we note that 
material properties of this form cannot be constructed exactly since the shear 
modulus /i0 becomes unbounded as r — )■ i?i (the inner boundary of the cloak) 
whilst the density tends to zero in this limit. 

In figure m we show an example of antiplane cloaking with a cloaking region 
defined by inner and outer radii i?i = 27r, R2 — 47r and the line source located at 
i?o = Stt, 80 = 0. This figure is plotted by solving (|2.3l) - (|2.5p which is achieved 
by noting that the solution to (|2.ip is (see lBj) 

W{R, 9) = W,{R, 6) = -^Ho(if 5) (2.6) 

4i/i 

where S = ^{X - X^f + [Y - Yaf and where Hq (isTi?) = H[,^^(Xi?) = io{KR) + 
iYo{KR) is the Hankel function of the first kind and order zero, noting that 
Jo and Yq are Bessel functions of the first and second kind respectively, also of 
order zero. We note that this solution retains the same form in the mapped con- 
figuration for R > R2 since the mapping in this region is the identity mapping. 
Therefore for r > R2 we have 

w{r,0) = ^IloiKs) (2.7) 

4:1 fj, 



where s = y/{x — Xq)^ + {y — Yq)"^, noting that the source location remains at 
X = Xo = RocosQo,y = Yq = i?osin6o- For < R < R2, corresponding to 
i?i < r < i?2, we use Graf's addition theorem with Rq > R2 (see (jB.Gp ) followed 
by the mapping p.2|) so that in the region Ri < r < R2 



C °° 



w{r,0) = -- V H„(i^i?o)J«(ifi?)e"'(®-®°) (2.8) 

4i/x ^-^ 



n— — 00 



where here i? = i?2(r - i?i)/(i?2 - Ri) and 6 = 6*. 

As discussed in the introduction, cloaking thus requires very special materi- 
als having properties such as (j2.5p above and the construction of such materials 
is non-trivial. In what follows we shall show that a similar cloaking effect can 
be obtained as that described above but with naturally obtained anisotropic, 
inhomogeneous shear moduli and a scalar, homogeneous density, thus not re- 
quiring the construction of special materials. This effect is due to nonlinear 
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Figure 1: A region of (nondimensionalized) radius Ka — 2tt is cloaked by the 
use of an elastic cloak in 2tt < Kr < Att whose properties (12.51) are obtained by 
applying the classical linear transformation (I2.2p to Helmholtz equation in two 
dimensions. Here the line source is located at Kr — KRq = Stt and Gq = and 
is shown as a small finite region to avoid the large values near its origin. 



elastic pre-stress. The corresponding equation still has the form of (|2.4p but 
now with 



fir{r) = A* 



M 



fJ-eir) =M 



M 



(2.9) 



where M = A^ — a^ relates the undeformed and deformed radii of the cylindrical 
cavities A and a respectively before and after pre-stress. We note further that 
the density remains uniform. 

3 Pre-stress 



Let us consider an isotropic incompressible neo-Hookean material of infinite 
extent with a cylindrical cavity of initial radius A at the origin. The constitutive 
behaviour of a neo-Hookean material is described by the strain energy function 
(Ogden 1997) 



M 



A^, 



W=^(/i-3) = ^(Af + A^ + A^-3) 



(3.1) 



where Ij and Xj,j = 1,2,3 are the principal strain invariants and principal 
stretches respectively, of the deformation which is to ensue. Stresses associated 



with the deformation are obtained via derivatives of this strain energy function, 
and are defined in p.lOp befow. We consider the initial deformation as depicted 
on the left of figure [2] which arises due to an increase in internal pressure inside 
the cavity region, leading to a cavity with larger radius. In order to describe 
this deformation mathematically we therefore write 

R^R{r), e = e, Z = z, (3.2) 

where {R, 0, Z) and (r, 6, z) are cylindrical polar coordinates in the undeformed 
and deformed configurations, denoted B and b respectively. The function R{r) 
will be determined from the radial equation of equilibrium and incompressibil- 
ity condition below. Note the convention introduced in p.2p . i.e. that upper 
case variables correspond to the undeformed configuration whilst lower case cor- 
responds to the deformed configuration. This notation is therefore analogous 
to that used for the untransformed and transformed coordinates for classical 
cloaking as described in section [2l Since we are interested in incremental per- 
turbations from the pre-stressed state, it will be convenient for us to derive 
equations in terms of coordinates in the deformed configuration. Position vec- 
tors in the undeformed (upper case) and deformed (lower case) configurations 
are respectively 

R{r) cos( 

i?(r) sin 61 | , x= ( rsin6i | . (3.3) 

Using p.2p . it can be shown that the principal stretches for this deformation 
in the radial, azimuthal and longitudinal directions, respectively, are 

dv 1 T 

The deformation gradient tensor F is defined by 

F = Gradx, (3.5) 

where Grad represents the gradient operator in X and here we have 

(i?/(r))-i 

r/R{r) | . (3.6) 

1 

For an incompressible material, we must have J — detF = 1, and so 






The differential equation p.7p can easily be solved to yield 



R{r) = V?^ + M, (3.8) 
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Figure 2: In an unbounded medium, a time-harmonic line source, located at 
(i?Oj0o) is indicated by crossed lines and we show schematically the direction 
of outward propagating waves and wave-fronts. This wave is scattered from a 
small cavity of radius A <^ A where A is the wavelength of the source, with 
internal pressure pm = in the undeformed configuration. When the pressure 
is increased, pin > 0, the cavity increases its radius to a such that a/A = 0(1). 
A region of inhomogeneous deformation develops around the deformed cavity. 



where Af is a constant, the form of which is defined by the corresponding radii 
of cavity boundaries, i.e. 

M = A'^-a'^, (3.9) 

with the equation linking a, A and internal pressure Pm given shortly. 

From Ogden (1997, 2007), the Cauchy and nominal stress tensors for an 
incompressible material are respectively given by 

dW dW 

where W is the neo-Hookean strain energy function introduced in p.ip . I is 
the identity tensor and Q is the scalar Lagrange multiplier associated with 
the incompressibility constraint. Note that the derivative with respect to F is 
defined component-wise by 



QF- 



(3.10) 



I 5F /„. OF,, 



(3.11) 



which is the convention introduced in Ogden (1997). The tensors F and T are 
diagonal, the non-zero components of the latter being given by 

Trr^^Ji\—^^\+Q. ^^6e = ^^[-^^J-^+Q, T,z = fi + Q. (3.12) 

The static equations of equilibrium are 

divT = 0, (3.13) 

where div signifies the divergence operator in x. These reduce to 

-Tf^ + -{Trr~Teg)=0, ^=0, ^=0. 3.14 

or r oO oz 

The second and third of these simply yield Q — Q{r). Integrating the first 
of p.l4p . using Trr and Tgg from p.l2p and applying the traction condition 
t\r=a = Tn|r=a = Pin^r where e^ is the unit vector in the radial direction 
and n = — e,. is the (outward) unit normal from the elastic material, so that 
Trr\r=a = ^Pin, wc find that thc radial stress is determined by 

r„(.-, +.,„ ^ I (m (1-1), log (Il±il) - log (:!)) . (3.15) 

This allows Q to be determined from (|3.12l) . 

Next using the condition that Trr — J- as r — >■ cxo in p.lSp yields the rela- 
tionship 

This nonlinear equation determines the deformed to undeformed radius ratio 
a/ A as a function of the scaled pressure Pin/M and we plot this relationship in 
figure [21 noting that moderate values oi pin/n can lead to relatively large values 
of a/ A. We note that for rubbery materials, the shear modulus is usually low, 
fi = O(105) Pa. 

4 Incremental deformations 

The antiplane line source will generate outgoing waves and subsequent scattered 
waves from the cavity. We now consider how to obtain the (incremental) equa- 
tion which governs wave propagation in the pre-stressed medium. This equation 
is derived by appealing to the theory of "small-on-large" (Ogden 1997, 2007). 
To this end let us consider incremental deformations from the (statically) de- 
formed body b. We consider a finite deformation of the original body i? to a 
new deformed state b which is close to the configuration b. The position vector 
in the new deformed state b is defined by x and we define 

u = x-x, |u| < 1 (4.1) 




Vin/^^ 



Figure 3: Plot of the deformed to undeformed radius ratio a/ A as a function of 
Pin/M- 



as the difference between position vectors in b and h. The term u is called an in- 
cremental displacement since b is close to 6. Since our line source is of antiplane 
type and is time-harmonic, we therefore assume that the total incremental de- 
formation is of the form 



u = ^[w{x, y)ez exp(— za;i)] , 



(4.2) 



where e^ is a unit basis vector in the z direction. On employing the theory of 
small-on-large, the technical details of which can be found in \^ we find that 
the governing equation of motion for w is 



19/ dw 

r dr \ ^ dr 



J.2 QQ2 



where 



/^rW = M 



M 



C 

+ puj'^w ^ —5{r - rQ)S(e - Bq) 



fJ-eir) = fJ. 



M 



(4.3) 



(4.4) 



and C is the force per unit length of the line source (as associated with the 
original line source). We note that as r — > cxo, (|4.3p reduces to 



1 d 



r dr \ dr 



dw 



\^+K'w^—Sir~ro)Si9~eo) (4.5) 



where K^ ~ puP' j [i is the wavenumber associated with the undeformed material. 
Note that as r — ^ cx) the deformation does not affect the wave and the equation 
is simply the scalar wave equation in the far field. 
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5 Cloaking via pre-stress 

5.1 Explicit solution of the incremental equation 

In [B] we derive the standard form of the scattered antiplane elastic wave field 
(|B.7|) - (|B.8|) which arises due to scattering from a cylindrical cavity (filled with 
inviscid fluid) in an unstressed medium, with a time-harmonic line source located 
at {Rq, 0o). In the undeformed configuration B we take KA <C 1 and therefore 
the leading order scattering coefficient is 0{{KA)'^) (Parnell and Abrahams 
2010), i.e. the scattered field is weak compared with the 0(1) incident field. 
When KA = 0{1) this scattered field is no longer weak in general. 

Let us now consider the pre-stressed problem in the deformed configuration 
b. The line source is thus moved to the location (rp, ^o) = (\/i?o — M, 6o) (this 
movement will itself be negligible provided it is far enough away from the cavity) 
and we retain the strength C of the line source associated with the undeformed 
configuration. Let us consider pressures Pm that are large enough in order to 
infiate the cavity to a radius a such that Ka = 0{1). 

The equation governing the incremental waves in the configuration b is (|4.3p . 
or explicitly 

(5.1) 

noting again that K^ = puP' j [i is the wavenumber associated with the unde- 
formed material. In this specific case of a neo-Hookean elastic material we can, 
in fact, solve this problem analytically and we show, perhaps surprisingly, that 
the scattering coefficients are completely unaffected by the pre-stress and there- 
fore even though the cavity has increased significantly in size, its scattering 
effect on incident waves remains weak. 
In (|5.ip . we introduce the mapping 

i?2 = r^ + M, e = 61 (5.2) 

and we note that this mapping corresponds exactly to the initial finite defor- 
mation (133) and dMl). Upon defining W(R, 6) = w{r(R), 0(0)), we find 

It is straightforward to show that 

-S{r-ro) = ^d{R-Ro) (5.4) 

ro Rq 

and hence (|5.3p becomes 

V'^W + K'^W =-^5{R-Ro)5{Q-Qo). (5.5) 

^ Ro 
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Wc conclude that the solution of (I5.5P is entirely equivalent to the solution of 
(|B.2[) obtained in[Bj The incident field is therefore given by (|B.5[) and (IB.6p . 
which thus gives rise to the scattered field associated with (j5.5p as (cf. (jB-Tp ) 



oo 

E 

n—~oo 



W,iR) = y (-i)"a«H„(i^ii')e"(«5-Q«). (5.6) 



Referring to |B1 and imposing the necessary traction free boundary condition on 
R — A, the scattering coefficients a„ are given by 

We now map back to the deformed configuration, using (|5.2I) in (|B.6p and 
(|5.6p . The incident and scattered fields in the deformed configuration arc there- 
fore 



(5.8) 
and 

oo 

«'-W= E (-')"a«H„ (i^v/r2 + Af) e™('^-''°) (5.9) 

n— — oo 

where we used Graf's addition theorem to obtain the form (j5.8p (Martin 2006). 

Perhaps most important to note is that the scattering coefficients a„ de- 
pend on the initial distance Rq between the centre of the cavity and the source 
location, and the undeformed cavity radius A and they are unaffected by the 
pre-stress, an — An- The wave field close to the cavity is affected by the pre- 
stress and in fact waves are bent around the object, whereas in the far-field 
the wave that is seen will be identical to the field that would be seen in that 
region when scattering takes place from the undeformed cavity. Since scattering 
is weak in the latter case {KA ^ 1), it will also be weak from the deformed 
(large) cavity {Ka = 0(1)). 

The analytical solution to the scattering problem above is used in Parnell 
and Abrahams (2011) to consider antiplane wave scattering in a pre-stressed 
microvoided composite with a neo-Hookean host phase. 

5.2 Cloaking 

In order to illustrate the above theory, consider the case where a/ A = 20 with 
KA = 2Tr/20, Ka = 27r. We set KRq = Btt, Oq = 0. As in figured on 
introducing the wavelength A via A = 2Tr/K, this choice ensures that A/ A = 20 
(weak scattering) whereas A/a = 1. Having the ratio A/a = 1 in an unstressed 
medium would normally result in strong scattering and the notable presence of 
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a scattered field and shadow region. Using pre-stress however, the scattering 
effect of a cylindrical cavity of radius a such that Ka — 0{1) is equivalent to 
that of the undeformed cavity with KA ^ 1, i.e. very weak. Note that an inner 
pressure of around Pm/A* = 3.5 ensures this deformation. 

We illustrate this situation with three plots of the total wave field in figure 
H On the left we show the total field W = Wi + Ws (defined in (|R6)) and (|BJ)) ) 
resulting from scattering from the small cylindrical cavity {KA = 27r/20) in an 
undeformed medium. Since A ^ ^, scattering is very weak and the scattering 
effect cannot be seen on the plot. In the middle plot we show the equivalent 
problem but with the pre-stressed cavity inflated from ^ to a so that Ka — 2n 
and the total displacement field plotted is w = Wi-\- Ws as defined in (|5.8p and 
(|5.9p . As shown theoretically above, the scattered field remains the same as for 
the undeformed problem and therefore even though Ka = 27r = 0{1), scattering 
is weak. This is in contrast to the plot on the right of the figure where we show 
the total field wave W = Wi + Ws for an undeformed cavity with KA = 2tt, 
illustrating the total wave field that would result in an unstressed medium in 
this scenario and the fact that it gives rise to strong scattering and a noticable 
shadow region. 
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Figure 4: Total antiplane wave-field due to a line source located at _Ro = 
Stt, 00 = in an undeformed configuration. The source is shown as a small 
finite region to avoid the large values near its origin. We show the total field 
in the case of a small cavity with KA = 27r/20 (left), a deformed cavity in the 
pre-stressed medium {A — >■ a) with Ka = 2tt (middle) and a large cavity with 
KA = 2-K in an unstressed medium (right). This illustrates the possibility of 
cloaking via pre-stress. 



6 Conclusions 

In this article we have shown how cloaking of objects from antiplane elastic 
waves can be achieved by nonlinear elastic pre-stress of a neo-Hookean material. 
In this case the deformation leaves the scattering coefficients unchanged and 
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thus provided the initial cavity is a weak scatterer, the deformed (large) cavity 
will also be a weak scatterer. In particular we note that the inhomogeneous 
anisotropic shear modulus that is found in the incremental equations arises 
naturally by virtue of the pre-stress. This means that no special metamaterials 
are required. Under this deformation the density also remains isotropic and 
homogeneous. That this was possible was due to the special transformation 
property of the incremental equations for neo-Hookean materials, in that the 
mapping which takes them back to the Helmholtz equation is exactly the same 
mapping (physical deformation) associated with the pre-stress. 

This same property does not hold for other nonlinear elastic materials (Mooney- 
Rivlin, Gent, Fung, etc.). This does not necessarily mean that a similar effect 
could not be achieved with these other materials however. If it can be shown 
that the scattering coefficients in the deformed configuration are changed only 
slightly from the scattering coefficients associated with the undeformed (weak 
scattering) cavity then similar effects would result. One should also consider 
finite regions to be pre-stressed (to form finite cloaks) which would require com- 
pressibility. 

Additionally, it appears that the notion of pre-stress permits active control: 
the size of the obstacle which can be cloaked is limited by the size of the pre- 
stressed (inflated) cavity. This is modified by changing the magnitude of the 
pressure inside the cavity. The continuous change in inhomogeneity induced 
by pre-stress also means that dispersive effects are not present, as they would 
be by the use of inhomogeneous metamaterials for example. This property is 
beneficial for broadband cloaking. 

We should note that there can be various problems associated with the 
stability of infiating cavities in elastic media. This issue may therefore require 
some consideration as to the construction of an appropriate elastic cloaking 
medium and in particular it may necessitate the use of elastomeric composites 
which can permit the inflation with greater ease. 

Finally, the obvious question is can this idea be generalized to other elas- 
todynamic problems? Cloaking in higher dimensional elastodynamic problems 
has proved difficult since Navier's equations are not invariant under coordinate 
transformations (Milton et al. 2006). Note that a two dimensional transforma- 
tion for in-plane elastic waves was considered in Brun et al. (2009) which led to 
invariance of the governing equations except that the transformed elastic modu- 
lus tensor did not possess the minor symmetries. Construction of materials with 
this property is non-trivial and motivated Norris and Shuvalov (2011) to consider 
the use of Cosserat materials which also do not possess these minor symmetries. 
Here, however we note that we have not used a coordinate transformation in 
order to motivate a cloak with inhomogeneous material properties - the map- 
ping which generates such inhomogeneous properies is a physical deformation 
associated with pre-stress. This naturally generates a transformed elastic mod- 
ulus tensor (the instantaneous modulus tensor M as defined in (jA.Bp ) which 
does not possess the minor symmetries. Therefore, provided we can ensure that 
the scattering coefficients in the pre-stressed configuration remain weak, even 
when Ka = 0(1) as was the case above, we can achieve approximate cloaking. 
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This clearly motivates studies on more general elastic waves in two and three 
dimensions. One could also hypothesize that pre-stress of nonlinearly elastic 
metamaterials could be useful in the cloaking context. 

A Summary of small-on-large theory 

For the application of small-on-large theory we follow the classical approach 
as in Ogden (1997, 2007) with slight modifications to notation. We consider 
incremental deformations from the deformed body b. To this end consider a 
finite deformation of the original body B to a new deformed state b which is 
close to the configuration b. Position vectors in the deformed states b and b are 
defined by x and x and we define 

u = x-x, |u| < 1 (A.l) 

as the difference between position vectors in b and b. Since b is close to b, u is 
called an incremental displacement which we assume is time-harmonic and of 
the antiplane type, i.e. 

u — 5R[w(x, y)ez exp{—iujt)] (A. 2) 

where ez is a unit vector in the 2;-direction in the body b, generated by a line 
source of amplitude C in the deformed configuration. 

Given the deformation gradient (13. 5|) , we can define the additional deforma- 
tion gradients f = gradx and F = Gradx = fF where grad and Grad denote 
the operation of taking the gradient in x and X respectively. Therefore 

r = Grad u = Grad(x - x) = (f - I)F (A.3) 

and similarly 7 = f — I, so that f = I + 7 and F = 7F. 

Furthermore, given the stresses (|3.10p in the configuration b, we can define 
the corresponding stresses T = T + r and S ~ S + s in 6 where t and s 
are the incremental Cauchy and nominal stresses respectively. Upon taking 
a Taylor expansion about the deformation state F, one can show that for an 
incompressible material 

s = L:-fF + qF-^-QF-^j (A.4) 

where q is the increment of Q and L — d^W/dF^ is defined componentwise via 
Lijki = d^yV/dFjidFik- The colon notation indicates A : b = Aijkihk- In (jA.4| 
we have neglected higher order terms, i.e. those of 0(7^). Defining the push- 
forward of the incremental nominal stress as ^ = Fs, from (|A.4I) this therefore 
has the form 

C = M : 7 + gl - Q7 (A.5) 

where the components of the instantaneous modulus tensor M are defined by 

^ijkl = r,^ a^ F.,„iFkn. (A.6) 
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Convenient forms for this tensor have been given by Ogden (1997) as 

Miijj = XiXjWij = Mjju, no sum on i, j, (A. 7) 

M,,,, = V AJ-A^ J ^- * ^ ■^' ^' ^ ^^' (A.8) 

[ XpV.i ~ AiAj Wy + A,W., i 7^ j, Ai = Aj, 

Mijj, = Mjiy =Mijy -AiWi, i^^j, (A.9) 

where Wi = dW /d\i and Wy = d^W/dXidXj. Analysis of the Cauchy stress 
gives rise to the relationship r = ^ + 7T for an incompressible material. 
The equations of motion in body b in terms of the Cauchy stress are 

d^V d^u ,, , 

where div denotes divergence in x, i.e. with respect to the configuration b. 
Equation (IA.10[) has this form since the body is incompressible {p = p) and 
tj = U + u where U is independent of time. Using the fact that divT = 0, it 
can be shown that for an incompressible medium 

dwt = div(r - 7T) = div C (A.ll) 

where we have used the relationship between ^ and t above. The incremental 
equations are therefore 

divC = P^ (A.12) 

together with appropriate forcing terms as required. Indeed, time-harmonic line 
source forcing in the z direction gives rise to incremental horizontally polarized 
shear waves of the form (jA.2[) . Two of the equations (jA.12|) imply that q = 
and the remaining single governing equation of motion is 

-I- K-) + -% + P^'y^ = -^(r - ro)6ie - 9o). (A.13) 

r or r oO tq 

where 

C,rz = Mrzrzlzr, Coz = MszBzJzB- (A. 14) 

In this case t^z = Crz and Tgz = Cez and so the instantaneous moduli M^zrz and 
MgzBz can be considered as pressure dependent anisotropic shear moduli which 
for an incompressible neo-Hookean material are 



Pr = Mrzrz = pX.,. = — 5 , pe ^ Mg^g^ ^ pXg - 



L\ r^ y ' r. .... r a iy^2^^j 

(A.15) 

Using (|A.14I) with (JA.ISP in (JA.13P gives the governing equation of motion in 
terms of the displacement w which is 

13/ , -,dw\ pg(r) d'^w , C* „ ^c/^ ^ -, / a -.^n 
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B Scattering of an antiplane elastic line source 
by a cylindrical cavity 

Consider an isolated cylindrical cavity of radius A, inside an isotropic homo- 
geneous elastic host region of infinite extent in all directions. We specify a 
Cartesian coordinate system {X,Y,Z), with origin at the centre of the cavity 
and whose Z axis runs parallel with the axis of the cylindrical cavity. Consider 
the problem of scattering by a time-harmonic line source (of small-amplitude) 
which is polarized in the Z direction and located at X = Xq, Y = Yo {a distance 
Rq away from the centre of the cylinder and at an angle 80 € [0, 2tt) subtended 
from the X-axis). This forcing therefore gives rise to linear elastic waves po- 
larized in the Z direction propagating in the XY plane, i.e. SH waves. The 
displacement field is therefore given by U = di.[W{X,Y)ez exp{—iujt)], where 
ez is a unit basis vector in the Z direction and lu is the circular frequency and 
where W{X, Y) is governed by (see e.g. Graff (1975)) 



,^)W - CS{X - Xo)S{Y - Yo) = ^S{R - Ro)S{Q - Go). (B.l) 



(AiV^+pcj^) 

Here we have defined the planar cylindrical polar coordinates R, O via the stan- 
dard relations X — R cos Q,Y — RcosQ and we note that Rq and 80 are thus 
defined by Xq = i?ocos8o,lo = i?ocos8o. We also point out that C is the 
force per unit length (in the Z direction) of the imposed line source and p and /i 
are the mass density and shear modulus of the host medium. We re- write (jB.ip 
in the form 

{V^+K^)W^-^S{R-Ro)S{e-eo), (B.2) 

where the wavenumber K is defined by K^ = puj'^ / p,. The traction- free boundary 
condition on the surface of the cavity at i? = ^ is 

dW 

We can write the solution of this problem in the form 

W = W, + Ws (B.4) 

where Wi and Ws represent incident and scattered (outgoing) fields from the 
cylindrical cavity. In particular, the incident field represents the solution of the 
inhomogeneous problem (jB.2p which is outgoing from the source location, i.e. 

W, = ^RoiKS) (B.5) 

where S = y/ {X - Xq)"^ + {Y - Yq)^ . We have defined Hq (-firs') = lii^\KS) = 
Jo{KS) + iYo{KS), the Hankel function of the first kind, noting that Jo and 
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Yq are Bessel functions of the first and second kind respectively, of order zero. 
Together with the time dependence in the problem, this ensures an outgoing 
field from the source. Graf's addition theorem allows us to write this field 
relative to the coordinate system {R, Q) centred at the origin of the cavity in 
the form (Martin 2006) 

W- = ^H (K<^^ = -^ X j^n=^oo^n{KRo)UKR)e^-(^-^"\ R < R,, 

(B.6) 

where H„ and J^ are respectively Hankel and Bessel functions of the first kind, 
and of order n. Upon writing the scattered field in the form 

oo 

Ws= Y. (-i)"A.H„(ifi?)e™(Q-®«) (B.7) 



n— — oo 



and using this together with (jB.4p and the first of (jB.6p in (IB. 31) we find that 
the scattering coefficients due to the incident wave from the line source are 

Let us now consider the limit as i?o ^ oo in (|B.5|) and the first of (|B.6p . Using 
the fact that (p. 364 of Abramowitz & Stegun (1965)) 



H^iKRo) ^ (-^)«^/-A_e»(i^«o^V4) (gj) 



as i?o -> oo and upon setting 



C = 2i^i^/2TTKRoe'^''^*^^"°\ (B.IO) 

we find that in this limit, 

oo 

W, ~ ^ i"J„(ifi?)e*"(®-^— ' =e*^(^'=°"^-'^+^"'"^— ), (B.ll) 



n— — oo 



which is an incident plane wave, propagating at an angle of incidence Oinc = 
&Q ~ TT £ [— 7r,7r) to the X axis (since 9o G [0,27r)). The scattering coefficients 
appearing in (jB.7p for this forcing arc thus 

" " KiKAy ^ ' 

where the superscript (pw) indicates a plane-wave forcing. These are the stan- 
dard scattering coefficients associated with plane wave scattering from a circular 
cylindrical cavity (see e.g. Martin (2006), p. 123, equation (4.5)), noting the form 
of the scattered field in (JB.7P with Qq = Oinc + tt. 
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